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Abstract Direct and resolved photon interactions are shown to be intimately related through 
the factorization mechanism. It is argued that in theoretically consistent analysis of jet 
production in 7p and ep collisions the LO resolved 7 contribution must be considered together 
with the NLO direct 7 component. Recent data from HERA therefore do not provide a direct 



evidence for the former component, but should rather be interpreted as a manifestation of 
^ ! the 0(a 2 a 2 s ) term in ep interactions. 
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In recent months an increasing flow of fresh data from HERA at DESY has brought 
several interesting new results on the structure of photon. Among them, the interpretation 
0] of data on two-jet events an as evidence for the resolved 7 contribution to 7p interac- 



[ tions has attracted considerable attention. I consider this interpretation as premature. Due 
to close relation between direct and resolved 7 components, based on the general factoriza- 
tion arguments, there is in fact no principal difference between these two contributions. I 
shall argue that in theoretically consistent treatment of 7p interactions, the usual resolved 
7 component should always be considered simultaneously with the NLO direct one, which, 
in fact, leads to kinematically very similar final state configurations. For the case of quarks 
inside the photon, the separation of 7p interactions into "direct" and "resolved" compo- 
nents is theoretically ambiguous and doesn't allow a unique interpretation of observed jet 
correlations as an evidence for the "resolved" photon. 

Let me start with an important comment on the photon structure functions. To the 
order 0(aa s ) the evolution equations for quark (of flavour i and charge e$), gluon and 
photon distribution functions inside the real photon read 

^ff = 0. D lh (x,M)=Hl-x) (1) 
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where ki(x) = 3e~(x + (1 — x) ) is the LO QED branching function corresponding to the 
vertex in Fig. la and P^l(z) are the analogous LO QCD branching functions. In (0-0) 
the argument (not written out explicitly) of the strong coupling constant a s , the so-called 
factorization scale M, is the same as the scale of the various distribution functions. In the 
above equation the term "quarks" denotes both quarks and antiquarks and correspondingly 
the summation in (Q) runs from 1 to 2nf. 

In many papers, e.g. [3-7], dealing with the photon distribution functions one finds the 
claim that the resolved 7 contribution, described by Feynman diagrams like those of Fig. 2a, 
is of the order 0(aa s ), despite the presence of two strong interaction vertices. This claim 
is based on incorrect analysis of the behaviour of D q / 7 (x, M) in the limit a s — > 0, which 
has lead to the wrong conclusion that it behaves like 0(a/a s ), the a s in the denominator 
cancelling one power of a s from the vertices. To show this, let us investigate, in LO QCD, 
the behaviour of D q i~ in weak coupling region. In the simplest case of the generic nonsinglet 
quark distribution function, D NS (x, M), the appropriate evolution equation reads || 

dD m (x,M) _ a u ^ i a s f 1 dy n ^ tjr \T>(o)( x 



dlnM 7r n Jx y \V J 

where 

Switching off QCD already at this stage, i.e. setting a s = on the r.h.s. of (Q), we get 

D^ D (x,M) = -k NS (x)\n^- } (6) 

where Mq is an arbitrary positive constant (M is by definition positive). Eq. (|6|) is, in fact, 
what we get by integrating the pole part 

t^nsOzO— (7) 

of the corresponding Feynman diagrams in Fig.2b, in the interval 1 1 |g (Mq , M 2 ). 

The claim that D NS (x, M) = C(a/a s ) or, in general, D q / ry (x, M) = 0(a/a s ), is based on 
the existence of the so-called "asymptotic pointlike" (or "anomalous" ) solution to the full 
evolution equation (f|), which is explicitly calculable in perturbative QCD 

DZ(x,M) = ^-a NS (x), (8) 
a s {M) 

where cxns^) is given as a solution to the equation 

aNsW = _|_, NsW + 2/;^p(;»g) 0Ns(!/) , (9 ) 

and which indeed behaves like 0(a/a s )\ However, as we shall see the existence of this 
solution provides no justification for the claim that the full solution of (HI) behaves like 
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0(a/a s ). Let me first recall the important fact H that this asymptotic pointlike solution 
doesn't exist alone but is, together with the solution of the corresponding homogeneous 
equation (describing the "hadronic" part of the photon), embedded in the general solution 
of the evolution equation (|4]). Converting both these components into moments for easier 
handling, we get for it the expression M: 



D NS (n,M) 



4tt 



a s (M) 



a 



l-2P(°)(n)//3 ' 



M M o)J 



a NS (n) + 



cts{M) 
a s (M ) 



-i -2F(°)(n)//3o 



which has been obtained as the sum of the asymptotic pointlike solution 
solution to the corresponding homogeneous equation 



D^(n,M)=D^(n,M ) 



a s {M) 
a s (M ) 



-|-2P(°)(n)//3 



D m (n,M ) 
(10) 

and the general 



taking into account (0). Mq in the above equations denotes the scale at which the boundary 
condition on D N $(x, Mq) is specified. 

In the weak coupling limit, i.e. for Aqcd — > 0, the second term in ([II]) goes over into 
M-independent D-^s(x, Mq), i.e. just the initial condition, while the dominant first term 
behaves as 
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2P(°)(n)//? )ln(M/M ) 



a s (M) 



ln(M/A QCD ) 



a NS (nJ 



l-2PW(n)/(3 )-k NS (n)ln^- (12) 

' 71 Mq 



and is thus very close to the simple expression (|6]), expected solely on the ground of QED 
vertex in Fig.l, with no QCD effects at all. The only difference with respect to (§) rests in 
the numerical factor (1 — 2P(°' (n) / (5q) containing the term proportional to P'°'(n), which 
is a consequence of the LO QCD corrections (like those of Fig. lb) to the basic QED vertex 
of Fig. la. As the limit flT2] ) is, as noticed already in H, finite, there is no justification to 
write D N g(x,M) = 0(a/a s ). The same is of course true for the singlet quark distribution 
function and for -D 9 / 7 as well. The preceding considerations clearly show how important 
it is to consider the asymptotic pointlike solution always in conjunction with the hadronic 
component of the photon. Recall, as another example, that D^ 7 (x, M) in (||) diverges badly 
at x — 0. Also in this case the problem is cured || in the sum ([II]). 

The crucial aspect of the above procedure is obviously the way the weak coupling limit is 
constructed. Let us compare the three above mentioned expressions (|5|,|3|,|TO"D for Dns(x, M) 
in this limit, obtained by sending Aqcd — ► 0. They differ in two aspects: the place where 
the limitting procedure was carried out and the interplay with the hadronic component. In 
the case of (^|), a s was set to zero in the evolution equation itself, leading to a well-defined, 
finite result for all M, which, however, carries no trace of the QCD corrections, present for 
a s ytz 0. For the asymptotic pointlike solution (§), considered separately from the hadronic 
part (|Tl"l) , the limit Aqcd — * leads to the result D^ s (x,M) — > oo for any x, M. The 



hadronic component (O), on the other hand, approaches a finite limit, given, as expected, 



by the initial distribution function D^(x,M )\ This is another signal that D^(x,M) has 
little physical meaning of its own but makes sense only in combination with the hadronic 
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part of the photon. In other words, taking first the limit Aqcd —> and then adding the 
two components leads to ill-defined result. 

If, on the other hand, we first add in (|T0|) the two mentioned components and only then 
take the limit Aqcd — ► we get not only a well-behaved result (0), but this finite result 
contains (in the term proportional to P(°'(n)) also a clear trace of the presence of higher 
order QCD corrections. Note that this finite limit, coming from the first term on the r.h.s. 
of (|10|), is actually independent of the parametrization of the hadronic part D^(x, M). 

In summary, the photon distribution function D q / y is clearly of the order O(o) and 
consequently the resolved photon contribution to, for instance, jet production in 7p collisions 
of the order 0(aa 2 ) and not, as claimed in [3-7], of the order 0(aa s ). It should therefore 
be added to the NLO direct 7 contribution, which is of the same order. 

The close relation between the LO resolved and NLO direct 7 contributions is in fact 
central to the very idea of factorization of parallel singularities. It also turns out that 
the factorization of parallel singularities provides another clear evidence that D q / 1 (x 1 M) = 
0(a). To demonstrate this assertion, let me discuss, at the NLO, the factorization of parallel 
singularities into the structure functions of beam particles in three distinct cases: 

p + p — > 2 jets + anything (13) 

7 + p — > 2 jets + anything (14) 

7 + p — > W + anything (15) 

Consider the basic QED vertex of Fig.l. Assuming that it is the lower, quark, leg which enters 
the interaction vertex in Fig. 2b, we encounter a singularity of the form \jt. Its factorization 
then amounts to dividing the whole integration range into the lower part | t |e (t min ,M 2 ) 
{t min is essentially an infrared regulator), put into the quark distribution function D q / 7 (x, M) 
and used in the LO resolved 7 contribution of Fig. 2a, and the upper one, | t |e [M 2 ,t max ), 
retained in the NLO hard scattering cross-section of the direct 7 subprocess of Fig. 2b. It 
seems obvious that due to the arbitrariness in the choice of the boundary value M 2 between 
what is included in the "resolved" and what is left in "direct" 7 component, these components 
have to be of the same order and thus D q / 1 (x, M) of the order 0(a). 

The general idea of factorization of parallel singularities [|l(J implies the following struc- 
ture of the cross-section for the process flT3|), described by Feynman diagrams in Fig. 3 Q 

gt( pp -> 2 jets) = W dxx f dx 2 D a/p ( Xl , M)D b/p (x 2 , M)cr^ rd (S, x u x 2 , M, /i), (16) 

a,b J ° J ° 

where the hard scattering cross-section cr^ rd for the parton level subprocess a-\-b — > c-\-d-\-- ■ ■ 
admits perturbation expansion in powers of a s (/x) at the hard scattering scale generally 
different from the factorization scale M: 



a h J rd (S lXl ,x 2l M^)= a^)K^(S,x 1 ,x 2 ) + aM^ u (S,x 1 ,x 2 ,M,iJ,) , (17) 



LO/c ™ \ 1 „.3/..\ NLO/ 



1 In Figures 3-5 only examples of typical Feynman diagrams, relevant to the discussed point, are shown. 
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where S is the total pp CMS energy, the sum in ( |I6| ) runs over all parton pairs a, b in 
the incoming hadrons. The functions k lo , k nlo are assumed to contain all appropriate 8- 
functions defining the kinematics of the final state jets in (p~3f) . Note that k lo (S, X\, x 2 ) is 
finite, M-independent function of these variables. In the above equations the factorization 
scale M is a free parameter, which separates short distances (large 1 1\ \ , \ t 2 | of the incoming 
partons in Fig. 3) from large distances, i.e. small virtualities, when both incoming partons a, b 
are close to their mass shell. The formal invariance of the factorization procedure with respect 
to the choice of M is guaranteed by the cancellation mechanism which, in the above process, 
works in such a way that the M- dependence of the LO contribution to fllBf ), contained in 
the M-dependent distribution functions D a / p (x, M) and described by Feynman diagrams of 
Fig. 3a, is cancelled to the NLO by the explicit dependence on M of k n , corresponding 
to diagrams of Fig.3b. Q This cancellation mechanism works separately for any two-parton 
final state in the subprocess (??). 

In the case of the photoproduction of jets ( |T4| ) the cancellation mechanism is modified 
by the presence of the inhomogeneous term in the evolution equation for quark distribution 
function D q /^(x, M). Let us consider, as an example, the two gluon final state of Fig.4a 
and, moreover, fix the lower parton leg, i.e. the parton coming from the target proton, 
to be quark of a particular flavour. The NLO QCD corrections to this subprocess are 
described, in the case of real emissions, by diagrams of Fig.4b. They cancel that part of 
the M-dependence of the LO diagrams of Fig. 4a, which is induced by the homogeneous 
part of the evolution equation for D g / 7 (x, M), i.e. the one involving the convolutions with 
the quark or gluon distribution functions. They, however, don't cancel the dependence of 
D q / y (x, M) on M induced by the inhomogeneous term in (^). This additional dependence 
on M clearly needs a "direct" photon in the initial state and is therefore cancelled by the 
"direct" 7 contribution, corresponding to the diagram in Fig. 4c. The cross-section of the 
2^3 subprocess 7+q^q+g+g has a parallel singularity, arising from the configuration 
where the quark-antiquark pair originating from the incoming photon is parallel to it, which 
has to be subtracted. Instead of ([16]) we thus have 

a( 7 +p -> 2 jets) = £ T dx 2 D b/p (x 2 , M)a^ rd (S, x 2 , M, //) 
b Jo 

+ E / dx ± T dx 2 D ah ( Xl , M)D b/p (x 2 , M)a^ rd (S, Xl ,x 2 , M, /1) (18) 

a,b J ° J ° 

where the direct 7 hard scattering cross-section 



a^riS, x 2 , M, n) = [aa 8 (ti)K%(S, x i} x 2 ) + aa 2 s (/i)< fe LU (5, Xl ,x 2 , M, fi)\ (19) 

contains the M-independent LO contribution, corresponding to diagram in Fig.4d, and M- 
dependent NLO part, which cancels the rest of the M-dependence of the LO "resolved" 7 
contribution induced by the inhomogeneous term in (H)! Thus the inclusion of the NLO 



2 Similar procedure can be formulated for the case of final state collinear singularities. These are ei- 
ther put into the appropriate fragmentation functions, or cancelled, via KLN theorem, in the sum over all 
indistinguishable partonic final states. 
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direct 7 contribution is neccessary for a theoretically consistent description of jet production 
in 7p collisions to the order 0(a.aP s ). As far concerned, one has to be careful as 

this quantity depends on the choice of the NLO branching function P^\z). Moreover, we 
cannot just subtract from the integrand of the appropriate Feynman diagrams the singular 
part, corresponding to the 1/t term, as the M-dependence enters through the upper limit 
on the interval of | t \, where this subtraction is carried out. This question is discussed in 
more detail in [|T3|j . 

At the 0(aa 2 ) order the factorization scale dependence doesn't cancel, but it is clear 
that the M-dependent NLO direct 7 term has to be included along with the LO resolved 
7 one, as they are both of the same order. For instance, in a recent analysis ]16| of the 
hadronic properties of the photon the total cross-section of 7p interactions is written as a 
sum of three terms, the first two of them corresponding to the "hadronic" and asymptotic 
pointlike components of the resolved 7 contribution, while the third one is the LO direct 
7 contribution. However, as emphasized above, any time the LO resolved 7 contribution 
is taken into account, so must also be, for theoretical consistency, the NLO direct 7 one! 
Moreover, as in |1| the scale M of D^ s (x,M) is taken very small (~ 0.5 GeV), the NLO 



direct 7 contribution evaluated at this scale may be quite large and its final state practically 
indistinguishable from the resolved 7 one. Although including this NLO direct 7 component 
means adding another positive contribution and thus further complicating the situation in 
ITjfl , there is no justification for neglecting it. 

Finally, at the 0(aa s ) order, only the M-independent direct 7 component contributes, 
which has no parallel singularity once the minimal transverse momentum of the produced 
jets is specified. The qq pair originating from the photon in Fig.4d can not be parallel and 
still produce nonzero p t at the lower vertex. 

In the case of W boson production in 7p collisions (p~5| ) the cancellation mechanism starts 
earlier than in the preceding subcase. The lowest order diagrams corresponding to (|i~5"D are, 
for direct as well as resolved 7 components, sketched in Figs. 5a, c. Both of these contributions 
are of the order 0(aaw)- As the produced W is massive, the qq pair originating in Fig.5c 
from the incoming photon can, contrary to diagram in Fig.4d, be parallel. Consequently, 
the parallel singularity of the 0(aaw) direct 7 term has to be subtracted, leading to the in- 
dependent finite part. The subtracted, M-dependent, part is included in the M-dependent 
photon distribution function D q / y (x, M), entering the LO contribution to the resolved 7 
component of Fig. 5a. Higher order QCD corrections (Fig. 5b) work as before. In this case the 
mechanism of factorization thus starts to operate already at the order 0(aaw) and provides 
another evidence against the claim that D q / y (x, M) = 0{a/a s ). Used in this process the 
latter behaviour would imply for the resolved component a res {^ + p — > W) = 0(aaw/oc s )- 
Consequently, the LO direct and resolved contributions would be of different orders in a s 
and, worst of all, a res would diverge for a s — > 0. 

So far I have discussed only the real photon interactions. In ep collisions the photon 
exchanged in NC interactions (Z° exchange can be neglected in the region of Q 2 considered) 
is virtual and one has to be careful to treat properly the dependence of photon distribution 
functions Dy^x, M,Q 2 ), i—q, q, g on the virtuality Q 2 of the exchanged photon fill . In 
this note I restrict my attention to the region of low Q 2 , where the virtual photon behaves 
to a very good approximation like the real one. The spectrum of these photons inside the 
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incoming electron is given by the Weizsacker- Williams approximation 



/ ~n " (l + (l-x) 2 \, Qs(l-x) , N 

where Qo denotes the maximum virtuality of the exchanged 7 taken into account. The "opti- 
mal" choice of Q (in the sense that the cross-section of the 2— >3 subprocess e+q/g— s-e+q+g/q 
behaves as l/Q 2 up to Q 2 = Qq) depends on the transverse momentum of produced partons 
and is roughly linear function thereof. My considerations concern this region. 

Up to now all theoretical analyses of data on jet production in ep collisions have included, 
along with the resolved 7 contribution of Fig.4a, also (and only) the LO direct 7 one of Fig.4d. 
These analyses show that around pi ~ 25 GeV, da res /dp t = da dir /dp t . As emphasized above, 
the 0(a 2 a 2 ) resolved 7 component should, however, always be considered simultaneously 
with NLO direct 7, which is of the same order. Although there exists WM a comprehensive 



analysis of direct 7 contributions to jet production in ep collisions up to the NLO, there 
is so far no theoretically consistent treatment of full jet production to the order 0(a 2 al), 
including LO resolved (Fig.4a) together with LO+NLO direct 7 contributions (Fig.4d,c). 
The results of [13] show that the ratio r = (da NLO / dp t ) / (da LO / dpt) is of the order of 2 in 



the region of pt around 20 GeV and that r grows as pt decreases so that the inclusion of NLO 
direct 7 contribution shifts the crossing point p\ to smaller values. It would be very useful 



to have the calculations of \TZ, I3| available down to p t around 10 GeV, where most of the 
available HERA data are, to see how much the situation will actually change by including 
the NLO direct 7 term. 

Moreover, in these circumstances the final state in the NLO direct 7 channel will often 
lead to kinematically very close final state configurations as the resolved one. This is due 
to the fact that the on-mass shell quark/ ant iquark, originating from the incoming photon 
and not participating in further interaction with the constituents of the proton, may fly 
close to the direction of the parent photon and thus is essentially distinguishable from the 
"remnant" jets of the resolved 7, Fig. 4a. For small p t the factorization mass M to be used in 
D q /y(x, M) should be roughly proportional to p t , which implies small opening angle between 
qq in the the direct 7 contribution of Fig.4c. For high p t the "remnant" jet from NLO direct 
component is expected to have bigger angle with respect to the beam (7) direction and thus 
to be distinguishable from the resolved 7 one. 

Let us now return to the interpretation of recent HERA data on two jet events. In 
only the LO direct and LO resolved 7 contributions were included. As theoretically consistent 
0(a 2 a 2 ) order analysis requires the inclusion of the NLO direct 7 term, the interpretation 
of the mentioned observation as an evidence for resolved 7 component is premature. What 
can be safely said is that the data require the presence of the 0(a 2 a 2 ) effects. Even this, 
however, is an important conclusion. 
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Figure captions 



Fig.l: Basic QED vertex (a) and example of lowest order QCD correction to it (b). In this 
and all following figures the big solid circles stand for parton distribution functions of the 
colliding hadrons (labelled by double lines) or the photon. 

Fig. 2: Examples of Feynman diagrams contributing to the resolved (a) and direct (b) 7 
photoproduction of jets. 

Fig. 3: Factorization for the two gluon final state in proton-proton collisions: the M-dependence 
of D q/p (x,M) used in C(a 2 s ) Fej r nman diagram of (a) is compensated by the explicit in- 
dependence of the 0(a^) hard scattering cross-section (b). 

Fig. 4: Factorization for the two gluon final state in 7p collisions: the M-dependence of 
D q / y (x, M) used in 0{aa 2 s ) Feynman diagram for resolved 7 (a) is compensated in part by 
the explicit M-dependence of the 0{aa z s ) resolved 7 hard scattering cross-sections (b) and 
in part by the M-dependence of the (D(aa 2 ) direct 7 hard scattering cross-section (c). In 
(d) the LO, 0(aa s ), direct 7 contribution is plotted for comparison. 

Fig. 5: Factorization in photoproduction of W boson: the M-dependence of D q / 1 (x, M) used 
in 0(aaw) resolved 7 contribution (a) is cancelled in part by the explicit M-dependence of 
0(aaw®s) hard scattering cross-sections (b) and in part by the explicit M-dependence of 
O(aaw) hard scattering cross-section (c). Dashed lines denote W boson. 
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